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VIDEO STRUCTURING BY PROBABILISTIC MERGING OF VIDEO 

SEGMENTS 

FIELD OF THE INVENTION 

5 . The invention relates generally to the processing and browsing of 

video material, and in particular to accessing, organizing and manipulating 
information from home videos. 

BACKGROUND OF THE INVENTION 

10 Among all the sources of video content, unstructured consumer 

video probably constitutes the content that most people are or would eventually be 
interested in dealing with. Organizing and editing personal memories by 
accessing and manipulating home videos represents a natural technological 
extension to the traditional still picture organization. However, although 

15 attractive with the advent of digital video, such efforts remain limited by the size 
of these visual archives, and by the lack of efficient tools for accessing, 
organizing, and manipulating home video information. The creation of such tools 
would also open doors to the organization of video events in albums, video baby 
books, editions of postcards with stills extracted from video data, multimedia 

2 0 family web-pages, etc. In fact, the variety of user interests suggests an interactive 

solution, which requires a minimum amount of user feedback to specify the 
desired tasks at the semantic level, and which provides automated algorithms for 
those tasks that are tedious or can be performed reliably. 

In commercial video, many moving image documents have story 
25 structures which are reflected in the visual content. In such situations, a complete 
moving image document is referred to as a video clip. The fundamental unit of 
the production of video is the shot, which captures continuous action. The 
identification of video shots is achieved by scene change detection schemes which 
give the start and end of each shot. A scene is usually composed of a small 

3 0 number of interrelated shots that are unified by location or dramatic incident. 

Feature films are typically composed of a number of scenes, which define a 
storyline for understanding the content of the moving image document 

In contrast with commercial video, unrestricted content and the 
absence of storyline are the main characteristics of home video. Consumer 
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contents are usually composed of a set of events, either isolated or related, each 
composed of one or a few shots, randomly spread along time. Such characteristics 
make consumer video unsuitable for video analysis approaches based on storyline 
models. However, there still exists a spatio-temporal structure, based on visual 
similarity and temporal adjacency between video segments (sets of shots) that 
appears evident after a statistical analysis of a large home video database. Such 
structure, essentially equivalent to the structure of consumer still images, points 
towards addressing home video structuring as a problem of clustering. The task at 
hand could be defined as the determination of the number of clusters present in a 
given video clip, and the design of an optimality criterion for assigning cluster 
labels to each frame/shot in the video sequence. This has indeed been the 
direction taken by most research in video analysis, even when dealing with 
storylined content. 

For example, in U.S. Patent No. 5,821,945, a technique is 
described for extracting a hierarchical decomposition of a complex video 
selection for browsing purposes, and combining visual and temporal information 
to capture the important relations within a scene and between scenes in a video. 
Thus, it is said, this allows the analysis of the underlying story structure with no a 
priori knowledge of the content. Such approaches perform video structuring in 
variations of a two-stage methodology: video shot boundary detection (shot 
segmentation), and shot clustering. The first stage is by far the most studied in 
video analysis (see, e.g., U. Gargi, R. Kasturi and S. H. Strayer, "Performance 
Characterization of Video-Shot-Change Detection Methods", IEEE CSVT, Vol. 
10, No. 1, February 2000, pp. 1-13). For the second stage, using shots as the 
fundamental unit of video structure, K-means, distribution-based clustering, and 
time-constrained merging techniques have all been disclosed in the prior art. 
Some of these methods usually require setting of a number of parameters, which 
are either application-dependent or empirically determined by user feedback. 

As understood in the prior art, hierarchical representations seem 
to be not only natural to represent unstructured content, but are probably the best 
way of providing useful non-linear interaction models for browsing and 
manipulation. Fortunately, as a byproduct, clustering allows for the generation of 
hierarchical representations for video content. Different models for hierarchical 
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organization have also been proposed in the prior art, including scene transition 
graphs (e.g., see the aforementioned U.S. Patent No. 5,821,945), and tables of 
contents based on trees, although the efficiency/usability of each specific model 
remains in general as an open issue. 

To date, only a few works have dealt with analysis of home video 
(e.g., see G. Iyengar and A. Lippman, "Content-based Browsing and Edition of 
Unstructured Video", IEEE ICME, New York City, August 2000; R. Lienhart, 
"Abstracting Home Video Automatically", ACM Multimedia Conference, 
Orlando, October, 1999, pp. 37-41; and Y. Rui andT. S. Huang, "A Unified 
Framework for Video Browsing and Retrieval", in A. C. Bovik, Ed., Handbook of 
Image and Video Processing , Academic Press, 1999). The work in the Lienhart 
article uses time-stamp information to perform clustering for generation of video 
summaries. Time-stamp information, however, might not always be available. 
Even though digital cameras include this information, users do not always use the 
time option. Therefore, a general solution cannot rely on this information. The 
work in the Rui and Huang article for generation of tables-of-contents, based on 
very simple statistical assumptions, was tested on some home videos with 
"storyline". However, the highly unstructured nature of home video makes the 
application of specific storyline models quite limited. With the exception of the 
Iyengar and Lippman article, none of the previous approaches have analyzed in 
detail the inherent statistics of such content. From this point of view, the present 
invention is more related to the work in N. Vasconcelos and A. Lippmann, "A 
Bayesian Video Modeling Framework for Shot Segmentation and Content 
Characterization", Proc. CVPR y 1997, that proposes a Bayesian formulation for 
shot boundary detection based on statistical models of shot duration, and to the 
work in the Iyengar and Lippmann article that addresses home video analysis 
using a different probabilistic formulation. 

Nonetheless, it is unclear from the prior art that a probabilistic 
methodology that uses video shots as the unit of organization could support the 
creation of a video hierarchy for interaction. In arriving at the present invention, 
statistical models of visual and temporal features in consumer video have been 
investigated for organization purposes. In particular, a Bayesian formulation 
seemed appealing to encode prior knowledge of the spatio-temporal structure of 



home video. In a departure from the prior art, the inventive approach described 
herein is based on an efficient probabilistic video segment merging algorithm 
which integrates inter-segment features of visual similarity, temporal adjacency, 
and duration in a joint model that allows for the generation of video clusters 
without empirical parameter determination. 

SUMMARY OF THE INVENTION 

The present invention is directed to overcoming one or more of the 
problems set forth above. Briefly summarized, according to one aspect of the 
present invention, a method for structuring video by probabilistic merging of 
video segments includes the steps of a) obtaining a plurality of frames of 
unstructured video; b) generating video segments from the unstructured video by 
detecting shot boundaries based on color dissimilarity between consecutive 
frames; c) extracting a feature set by processing pairs of segments for visual 
dissimilarity and their temporal relationship, thereby generating an inter-segment 
visual dissimilarity feature and an inter-segment temporal relationship feature; 
and d) merging video segments with a merging criterion that applies a 
probabilistic analysis to the feature set, thereby generating a merging sequence 
representing the video structure. In the preferred embodiment, the probabilistic 
analysis follows a Bayesian formulation and the merging sequence is represented 
in a hierarchical tree structure that includes a frame extracted from each segment. 

As described above, this invention employs methods for consumer 
video structuring based on probabilistic models. More specifically, the invention 
proposes a novel methodology to discover cluster structure in home videos, using 
video shots as the unit of organization. The methodology is based on two 
concepts: (i) the development of statistical models (e.g., learned joint mixture 
Gaussian models) to represent the distribution of inter-segment visual similarity 
and an inter-segment temporal relationship, including temporal adjacency and 
duration of home video segments, and (ii) the reformulation of hierarchical 
clustering (merging) as a sequential binary classification process. The models are 
used in (ii) in a probabilistic clustering algorithm, for which a Bayesian 
formulation is useful since these models can incorporate prior knowledge of the 
statistical structure of home video, and which offers the advantages of a 
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principled methodology. Such prior knowledge can be extracted from the detailed 
analysis of the cluster structure of a real home video database. 

The video structuring algorithm can be efficiently implemented 
according to the invention and does not need any ad-hoc parameter determination. 
5 As a byproduct, finding video clusters allows for the generation of hierarchical 
representations for video content, which provide nonlinear access for browsing 
and manipulation. 

A principal advantage of the invention is that, based on the 
performance of the methodology with respect to cluster detection and individual 
1 0 shot-cluster labeling, it is able to deal with unstructured video and video with 

unrestricted content, as would be found in consumer home video. Thus, it is the 
first step for building tools for a system for the interactive organization and 
retrieval of home video information. 

As a methodology for consumer video structuring based on a 
15 Bayesian video segment merging algorithm, another advantage is that the method 
automatically governs the merging process, without empirical parameter 
determination, and integrates visual and temporal segment dissimilarity features 
in a single model 

Furthermore, the representation of the merging sequence by a tree 
2 0 provides the basis for a user-interface that allows for hierarchical, non-linear 
access to the video content. 

These and other aspects, objects, features and advantages of the - 
present invention will be more clearly understood and appreciated from a review 
of the following detailed description of the preferred embodiments and appended 

2 5 claims, and by reference to the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram providing a functional overview of video 
structuring according to the present invention. 

3 0 FIG. 2 is a flow graph of the video segment merging stage shown 

in Figure 1 . 

FIG. 3 is a distribution plot of consumer video shot duration for a 
group of consumer images. 
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FIG. 4 is a scatter plot of labeled inter-segment feature vectors 
extracted from a home video. 

FIG. 5 is a tree representation of key frames from atypical home 

video. 

DETAILED DESCRIPTION OF THE INVENTION 

Because video processing systems employing shot detection and 
cluster analysis are well known, the present description will be directed in 
particular to attributes forming part of, or cooperating more directly with, a video 
structuring technique in accordance with the present invention. Attributes not 
specifically shown or described herein may be selected from those known in the 
art. In the following description, a preferred embodiment of the present invention 
would ordinarily be implemented as a software program, although those skilled in 
the art will readily recognize that the equivalent of such software may also be 
constructed in hardware. Given the system as described according to the invention 
in the following materials, software not specifically shown, suggested or 
described herein that is useful for implementation of the invention is conventional 
and within the ordinary skill in such arts. If the invention is implemented as a 
computer program, the program may be stored in conventional computer readable 
storage medium, which may comprise, for example; magnetic storage media such 
as a magnetic disk (such as a floppy disk or a hard drive) or magnetic tape; optical 
storage media such as an optical disc, optical tape, or machine readable bar code; 
solid state electronic storage devices such as random access memory (RAM), or 
read only memory (ROM); or any other physical device or medium employed to 
store a computer program. 

Accessing, organizing and manipulating personal memories stored 
in home videos constitutes a technical challenge, due to its unrestricted content, 
and the lack of clear storyline structure. In this invention, a methodology is 
provided for structuring of consumer video, based on the development of 
parametric statistical models of similarity and adjacency between shots, the unit 
of visual information in consumer video clips. A Bayesian formulation for 
merging of shots appears as a reasonable choice as these models can encode prior 
knowledge of the statistical structure of home video. Therefore, the methodology 



is based on shot boundary detection and Bayesian segment merging. Gaussian 
Mixture joint models of inter-segment visual similarity, temporal adjacency and 
segment duration -learned from home video training samples using the 
Expectation-Maximization (EM) algorithm- are used to represent the class- 
5 conditional densities of the observed features. Such models are then used in a 
merging algorithm consisting of a binary Bayes classifier, where the merging 
order is determined by a variation of Highest Confidence First (HCF), and the 
Maximum a Posteriori (MAP) criterion defines the merging criterion. The 
merging algorithm can be efficiently implemented by the use of a hierarchical 
q 10 queue, and does not need any empirical parameter determination. Finally, the 
^ representation of the merging sequence by a tree provides the basis for a user- 

rU interface that allows for hierarchical, non-linear access to the video content. 

« Referring first to Figure 1 , the video structuring method is shown 

^7 to operate on a sequence of video frames stage 8 obtained from an unstructured 

s 15 video source, typically displaying an unrestricted content, such as found in 

i i 

5j consumer home videos. The salient features of the video structuring method 

4f according to the invention can be concisely summarized in the following four 

q stages (which will be subsequently described in later sections in more detail): 

1) The Video Segmentation Stage 10: Shot detection is computed by 
20 adaptive thresholding of a histogram difference signal. 1-D color histograms are 

computed in RGB space, with N = 64 quantization levels for each band. The LI 
metric is used to represent the dissimilarity dc (t,t+l) between two consecutive 
frames. As a post-processing step, an in-place morphological hit-or-miss 
transform is applied to the binary signal with a pair of structuring elements that 

2 5 eliminate the presence of multiple adjacent shot boundaries. 

2) The Video Shot Feature Extraction Stage 12: It is known in the art that 
visual similarity is not enough to differentiate between two different video events 
(e.g., see the Rui and Huang article). Both visual similarity and temporal 
information have been used for shot clustering in the prior art. (However, the 

3 0 statistical properties of such variables have not been studied voider a Bayesian 

perspective.) In this invention, three main features in a video sequence are 
utilized as criteria for subsequent merging: 



• Visual similarity is described by the mean segment histogram 
that represents segment appearance. The mean histogram 
represents both the presence of the dominant colors and their 
persistence within the segment. 

• Temporal separation between segments is a strong indication of 
their belonging to the same cluster. 

• Combined temporal duration of two individual segments is also 
• a strong indicator about their belonging to the same cluster 

(e.g., two long shots are not likely to belong to the same video 
cluster). 

3) The Video Segment Merging Stage 14: This step is carried out by 
formulating a two-class (merge/not merge) pattern classifier based on Bayesian 
decision theory. Gaussian Mixture joint models of inter-segment visual 
similarity, temporal adjacency and segment duration -learned from home video 
training samples using the Expectation-Maximization (EM) algorithm- are used to 
represent the class -conditional densities of the observed features. Such models are 
then used in a merging algorithm comprising a binary Bayes classifier, where the 
merging order is determined by a variation of Highest Confidence First (HCF), 
and the Maximum a Posteriori (MAP) criterion defines the merging criterion. The 
merging algorithm can be efficiently implemented by the use of a hierarchical 
queue, and does not need any empirical parameter determination. A flow graph of 
the merging procedure is given in Figure 2 and will be described in further detail 
later in this description. 

4) The Video Segment Tree Construction Stage 16: The merging 
sequence, i.e. a list with the successive merging of pairs of video segments, is 
stored and used to generate a hierarchy, whose merging sequence is represented 
by a binary partition tree 1 8. Figure 5 shows a tree representation from a typical 
home video. 

1. An Overview of the Approach 

Assume a feature vector representation for video segments, i.e., 
suppose that a video clip has been divided into shots or segments (where a 
segment is composed of one or more shots), and that features that represent them 
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have been extracted. Any clustering procedure should specify mechanisms both 
to assign cluster labels to each segment in the home video clip and to determine 
the number of clusters (where a cluster may encompass one or more segments). 
The clustering process needs to include time as a constraint, as video events are of 
limited duration (e.g., see the Rui and Huang article). However, the definition of 
a generic generative model for intra-segment features in home videos is 
particularly difficult, given their unconstrained content Instead, according to the 
present invention, home video is analyzed using statistical inter-segment models. 
In other words, the invention proposes to build up models that describe the 
properties of visual and temporal features defined on pairs of segments. Inter- 
segment features naturally emerge in a merging framework, and integrate visual 
dissimilarity, duration, and temporal adjacency. A merging algorithm can be 
thought of as a classifier, which sequentially takes a pair of video segments and 
decides whether they should be merged or not. Let s t and Sj denote the i-th and 

j - th video segments in a video clip, and let s be a binary random variable (r.v.) 
that indicates whether such pair of segments correspond to the same cluster and 
should be merged or not. The formulation of the merging process as a sequential 
two-class (merge/not merge) pattern classification problem allows for the 
application of concepts from Bayesian decision theory (for a discussion of 
Bayesian decision theory, see, e.g., R.O. Duda, P.E. Hart and D.G. Stork, Pattern 
Classification , 2 nd ed, John Wiley and Sons, 2000). The Maximum a Posteriori 
(MAP) criterion establishes that given an n-dimensional realization x v of an r.v. x 

(representing inter-segment features and detailed later in the specification), the 
class that must be selected is the one that maximizes the a posteriori probability 
mass function of s given x , i.e., 



s* = argmax Pr(^ | x) 



s 



By Bayes rule, 



Pr(^ | x) = 



p(x\e)Pr(e) 
P(x) 



# . -10- • 



where p(x \ e) is the likelihood of x given s , and Pr(e) is the prior of e , and 
p(x) is the distribution of the features. The application of the MAP principle can 
then be expressed 




1, p(x\e = l)Pr(> = i)>p( x \e = 0)Pr(e = 0) 
0 otherwise 



or in standard hypothesis testing notation, the MAP principle can be expressed as 

p(x\e = \)Vt( < e = \) > p(x\e = 0)Viis = G) 
H 0 

where H x denotes the hypothesis that the pair of segments should be merged, and 
H 9 denotes the opposite. With this formulation, the classification of pairs of shots 
is performed sequentially, until a certain stop criteria is satisfied. Therefore, the 
tasks are the determination of a useful feature space, the selection of models for 
the distributions, and the specification of the merging algorithm. Each of these 
steps are described in the following sections of the description. 

2. Video Segmentation 

To generate the basic segments, shot boundary detection is 
computed in stage 1 0 by a series of methods to detect the cuts usually found in 
home video (see, e.g., U. Gargi, R. Kasturi and S. H. Strayer, "Performance 
Characterization of Video-Shot-Change Detection Methods", IEEE CSVT y Vol. 
10, No. 1, February 2000, pp. 1-13). Over-segmentation due to detection errors 
(e.g. due to illumination or noise artifacts) can be handled by the clustering 
algorithm. Additionally, videos of very poor quality are removed. 

In implementing a preferred embodiment of the invention, shot 
detection is determined by adaptive thresholding of a histogram difference signal. 
1-D color histograms are computed in the RGB space, with AT = 64 quantization 
levels for each band. Other color models (LAB or LUV) could be used, and might 
provide better shot detection performance, but at increased computational cost. 
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The LI metric is used to represent the color dissimilarity dc(t,t+l) between two 
consecutive frames: 

«*c('.'+i) = 2X-^ii 

*=i 

where h k t denotes the value of the k-th bin for the concatenated RGB histogram of 
frame t. The 1-D signal dc is then binarized by a threshold that is computed on a 
sliding window centered at time t of length fr/2, where fr denotes the frame rate. 




1 d c (t)> ju d (t) + ka d (t) 
0 otherwise 



where (t) denotes the mean of dissimilarities computed on the sliding window, 
cr d (t) denotes the mean absolute deviation of the dissimilarity within the 
window, which is known to be a more robust estimator of the variability of a data 
set around its mean, and A: is a factor that sets the confidence interval for 
determination of the threshold, set in the interval. Consecutive frames are 
therefore deemed to belong to the same shot if s(t) = 0, and a shot boundary 
between adjacent frames is identified when s(t)=l . 

As a post-processing step, an in-place morphological hit-or-miss 
transform is applied on the binary signal with a pair of structuring elements that 
eliminate the presence of multiple adjacent shot boundaries, 

b(t) = s(t)®(e l (t),e 2 (t)) 

where® denotes hit-or-miss, and the size of the structuring elements is based on 
the home video shot duration histograms (home video shots are unlikely to last 
less than a few seconds), and it is set to fr/2 (see Jean Serra: Image Analysis and 
Mathematical Morphology, Vol. 1, Academic Press, 1982). 

3. Video Inter-segment Feature Definition 
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A feature set for visual dissimilarity, temporal separation and 
accumulated segment duration is generated in the video shot feature extraction 
stage 12. Both visual dissimilarity and temporal information, particularly 
temporal separation, have been used for clustering in the past. In the case of 
visual dissimilarity, and in terms of discerning power of a visual feature, it is clear 
that a single frame is often insufficient to represent the content of a segment. 
From the several available solutions, the mean segment color histogram is 
selected to represent segment appearance, 

1 e ' 
m. = 

where h t denotes the t-th color histogram, and m. denotes the mean histogram of 
segment s n each consisting of M t . = e { —b. + 1 frames (b. and e i denote the 
beginning and ending frame of segment s i ). The mean histogram represents both 
the presence of the dominant colors and their persistence within the segment. The 
LI norm of the mean segment histogram difference is used to visually compare a 
pair of segments i and j, 

B 

where denotes visual dissimilarity between segments i and j 9 B is the number 
of histogram bins, m. k is the value of the k-th bin of the mean color histogram of 
segment s. , and mju is the value of the k-th bin of the mean color histogram of 
segment sj. 

In the case of temporal information, the temporal separation 
between segments s, and , which is a strong indication of their belonging to the 
same cluster, is defined as 

fi^nnnQei-bjUej-b, |)(1-^) 

where s v denotes a Kronecker's delta, b n e { denote first and last frames of 
segments , and bj,ej denote first and last frames of segment sj. 

Additionally, the accumulated segment (combined) duration of two 
individual segments is also a strong indication about their belonging to the same 
cluster. Fig. 3 shows the empirical distribution of home video shot duration for 
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approximately 660 shots from a database with ground-truth, and its fitting by a 
Gaussian mixture model (see next subsection). (In Figure 3, the empirical 
distribution, and an estimated Gaussian mixture model consisting of six 
components, are superimposed. Duration was normalized to the longest duration 
found in the database (580 sec.).) Even though videos correspond to different 
scenarios and were filmed by multiple people, a clear temporal pattern is present 
(see also the Vasconcelos and Lippmann article). The accumulated segment 
duration r ti is defined as 

r & = card (s ( ) + card Qij ) 

where card{s) denotes the number of frames in segment s . 
4. Modeling of Likelihoods and Priors 

The statistical modeling of the inter-segment feature set is 
generated in the video segment merging stage 14. The three described features 
become the components of the feature space x , with vectors x = (a,fi,r). To 
analyze the separability of the two classes, Fig. 4 shows a scattering plot of 4000 
labeled inter-segment feature vectors extracted from home video. (Half of the 
samples correspond to hypothesis H x (segment pair belongs together, labeled with 
light gray), and the other half to H 0 (segment pair does not belong together, 
labeled with dark gray). The features have been normalized. ) 

The plot indicates that the two classes are in general separated. A 
projection of this plot clearly illustrates the limits of relying on pure visual 
similarity. A parametric mixture model is adopted for each of the class- 
conditional densities of the observed inter-segment features, 

P (x i s,e) = 21 Pr ( c = 0p(* \sA ) 
t-i 

where K £ is the number of components in each mixture, Pr(c = /) denotes the prior 
probability of the / - th component, p(x \e > e i ) is the / - th pdf parameterized by 9, , and 
e = {Pr(c),{6> }} represents the set of all parameters. In this invention, we assume 
multivariate Gaussian forms for the components of the mixtures in d -dimensions 
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so that the parameters 0, are the means p. and covariance matrices r, (see Duda et 
al., Pattern Classification , op. cit). 

The well-known expectation-maximization (EM) algorithm constitutes the 
standard procedure for Maximum Likelihood estimation (ML) of the set of 
parameters © (see A. P. Dempster, N. M. Laird and D. B. Rubin, "Maximum 
Liklihood from Incomplete Data via the EM Algorithm", Journal of the Royal 
Statistical Society, Series B, 39:1-38, 1977). EM is a known technique for finding 
ML estimates for a broad range of problems where the observed data is in some 
sense incomplete. In the case of a Gaussian Mixture, the incomplete data are the 
unobserved mixture components, whose prior probabilities are the parameters 
(Pr(c)} . EM is based on increasing the conditional expectation of the log- 
likelihood of the complete data given the observed data by using an iterative hill- 
climbing procedure. Additionally, model selection, i.e., the number of 
components of each mixture can be automatically estimated using the Minimum 
Description Length (MDL) principle (see J. Rissanen, "Modeling by Shortest 
Data Description", Automatica, 14:465-471, 1978). 



increasing the conditional expectation of the log-likelihood of the complete data 
Y given the observed data X = {jc, x N } : 



by using an iterative hill-climbing procedure. In the previous equation, X - h(Y) 
denotes a known many-to-one function (for example, a subset operator), x 
represents a sequence or vector of data, and p is an superscript that denotes the 
iteration number. The EM algorithm iterates the following two steps until 
convergence to maximize Q{&) : 



The general EM algorithm, valid for any distribution, is based in 



Q{0 | e ip) ) = £{log p(Y | 0) | x,0 ip) } 
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E-step: Find the expected likelihood of the complete data as a function of 0 , 

Q(0\0 ip) ). 

M-step: Re-estimate parameters, according to 
g(P+D = a rgmax0(0|0 o ' ) ) 

In other words, firstly estimate values to fill in for the incomplete 
data in the E-Step (using the conditional expectation of the log-likelihood of the 
complete data given the observed data, instead of the log-likelihodd itself). Then, 
compute the maximum likelihood parameter estimate using in the M-step, and 
repeat until a suitable stopping criterion is reached. EM is an iterative algorithm 
that converges to a local maximum of the likelihood of the sample set. 

For the specific case of multivariate Gaussian models, the 
complete data is given by Y = (X, I) , where / indicates the Gaussian component 
that has been used in generating each sample of the observed data. Element-wise, 
y = (x,i), i e {1,..., K £ } . In this case, EM takes a further simplified form: 

E-step: For all N training samples, and for all mixture components, compute the 
probability that Gaussian i fits the sample x ; given the current estimation & ip) , 

M-step: Re-estimate parameters, 
^ =^f t p(i\x j ,e,®^) 



-16- 



N 




N 



^p(i\xj,eM p) ) 



N 



y Zp(!\x Jt s,ei»Xxj-/'F i >Kx J -n^Y 




N 



2>O|*,,*,0<*>) 



The mean vectors and covariance matrices for each of the mixture components 
must be initialized in the first place. In this implementation, the means are 
initialized using the traditional K-means algorithm, while the covariance matrices 
are initialized with the identity matrix. As other hill climbing methods, data- 
driven initialization usually performs better than pure random initialization. 
Additionally, on successive restarts of the EM iteration, a small amount of noise 
is added to each mean, to diminish the procedure to be trapped in local maxima. 

The convergence criterion is defined by the rate of increase on the 
log-likelihood of the observed data in successive iterations, 



N 



\o%L{®\X) = \ogY{p{xj\e,®) 



i.e., the EM iteration is terminated when 



logZ(Q (/?+1) | X) - log£(0 (p) | X) 
\ogL{@ ip) \X) 



<10 



-2 



The specific model, i.e., the number of components K s of each 
mixture is automatically estimated using the Minimum Description Length 
(MDL) principle, by choosing 
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Kl = arg max ( log L(0 | X) - ^Mog N) 
k s 2 

where L( ) denotes the likelihood of the training set, and n Kg is the number of 

parameters needed for the model, which for a Gaussian mixture is equal to 

■ x< -(*.-i) + *.rf + jr.««!Ltfi 

When two models fit the sample data in a similar way, the simpler model (smaller 
K s ) is chosen. 

Instead of imposing independence assumptions among the 
variables, the full joint class-conditional pdfs are estimated. The ML estimation of 
the parametric models for P (x \ € = 0) and p(x \ € = i) , by the procedure just 
described, produces probability densities represented by ten components in both 
cases, respectively. 

In the Bayesian approach, the prior probability mass function 
Pr(£-) encodes all the previous knowledge at hand about the specific problem. In 
this particular case, this represents the knowledge or belief about the merging 
process characteristics (home video clusters mostly consist of only a few shots). 
There exist a variety of solutions that can be explored: 

The simplest assumption is Pr(£ = 0) = Pr(^ = 1) = 1 / 2 , which turns the MAP 
criterion into the ML criterion. 

The priors themselves can be ML-estimated from training data (see Duda et 
al., Pattern Classification , op, cit). It is straightforward to show that, assuming 
that the TV are independent, the ML estimator of the priors is 

Pr(ff = e) = ^£i(e,*) 

where i(e,k) is equal to one if the k-th training sample belongs to the class 
represented by s = e , e e {0,1} , and zero otherwise. In other words, the priors are 
simply weights determined by the available evidence (the training data). 
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The dynamics involved in the merging algorithm (presented in the following 
section) also influences the prior knowledge in a sequential manner (it is 
expected that more segments will be merged at the beginning of the process, 
and less at the end). In other words, the prior can be dynamically updated 
based on this rationale. 

5. Video Segment Clustering 

The merging algorithm is implemented in the video segment 
merging stage 14. Any merging algorithm requires three elements: a feature 
model, a merging order, and a merging criterion (L. Garrido, P. Salembier, D. 
Garcia, "Extensive Operators in Partition Lattices for Image Sequence Analysis", 
Sign. Proa, 66(2): 157-180, 1998). The merging order determines which clusters 
should be probed for possible merging at each step of the process. The merging 
criterion decides whether the merging should occur or not. The feature model of 
each cluster should be updated if a merging occurs. The present video segment 
clustering method uses this general formulation, based on the statistical inter- 
segment models developed in the previous section. In the present algorithm, the 
class-conditionals are used to define both the merging order and the merging 
criterion. 

Merging algorithms can be efficiently implemented by the use of 
adjacency graphs and hierarchical queues, which allow for prioritized processing. 
Elements to be processed are assigned a priority, and introduced into the queue 
according to it. Then, the element that is extracted at each step is the one that has 
the highest priority. Hierarchical queues are now traditional tools in mathematical 
morphology. Their use in Bayesian image analysis first appeared in C. Chou and 
C. Brown, "The Theory and Practice of Bayesian Image Labeling", IJCV, 4, pp. 
185-210, 1990, with the Highest Confidence First (HCF) optimization method 
The concept is intuitively appealing: at each step, decisions should be made based 
on the piece of information that has the highest certainty. Recently, similar 
formulations have appeared in morphological processing. 

As shown in Figure 2, the segment merging method comprises two 
stages: a queue initialization stage 20 and a queue updating/depletion stage 30. 



The merging algorithm comprises a binary Bayes classifier, where the merging 
order is determined by a variation of Highest Confidence First (HCF), and the 
Maximum a Posteriori (MAP) criterion defines the merging criterion. 

Queue initialization. At the beginning (22) of the process, inter- 
shot features x {J are computed for all pairs of adjacent shots in the video. Each 
feature is introduced (24) in the queue with priority equal to the probability of 
merging the corresponding pair of shots, Pr(* = i\x v ) , 

Queue depletion/updating. The definition of priority allows 
making decisions always on the pair of segments of highest certainty. Until the 
queue is empty (32), the procedure is as follows: 

1 . In the element extraction stage 34, extract an element (pair of 
segments) from the queue. This element is the one that has the highest priority. 

2. Apply the MAP criterion (36) to merge the pair of segments, 

i.e., 

I e - l)Pr(* = 1) > p{x ¥ \e = 0)Pr<<r = 0) 

3. If the segments are merged (the path 38 indicating the 
application of hypothesis H x ), update the model of the merged segment in 
segment model updating stage 40, then update the queue in the queue updating 
stage 42 based on the new model, and go to step 1 . Otherwise, if the segments are 
not merged (the path 44 indicating the application of hypothesis H 0 ), go to step 1. 

When a pair of segments is merged, the model of the new segment 
s. is updated by 

m.. = (card (s. )m i + card (s^ ) ligard (s ( )+card £ 3 )) 

b. = min(6 ( ,^) 

e. = max(e.,e > ) 
card ) = card (y . ) -heard (Sj ) 

After having updated the model of the (new) merged segment, four 
functions need to be implemented to update the queue: 
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1 . Extraction from the queue of all those elements that involved 
the originally individual (now merged) segments. 

2. Computation of new inter-segment features x = (a,>f?,r)using the 
updated model. 

3. Computation of new priorities Pr(* = 1 1 x {J ) . 

4. Insertion in the queue of elements according to new priorities. 
Note that, unlike many previous methods (such as described in the 

Rui and Huang article), this formulation does not need any empirical parameter 
determination. 

The merging sequence, i.e., a list with the successive merging of 
pairs of video segments, is stored and used to generate a hierarchy. Furthermore, 
for visualization and manipulation, after emptying the hierarchical queue in the 
merging algorithm, further merging of video segments is allowed to build a 
complete merging sequence that converges into a single segment (the whole video 
clip). The merging sequence is then represented by a partition tree 18 (Figure 1), 
which is known to be an efficient structure for hierarchical representation of 
visual content, and provides the starting point for user interaction. 

6. Video Hierarchy Visualization. 

An example of a tree representation stage 50 appears in Fig. 5. 
A prototype of an interface to display the tree representation of the analyzed home 
video may be based on key frames, that is, a frame extracted from each segment. 
A set of functionalities that allow for manipulation (correction, augmentation, 
reorganization) of the automatically generated video clusters, along with cluster 
playback, and other VCR capabilities may be applied to the representation. The 
user may parse the video using this tree representation, retrieve preview clips and 
do video editing. 

Queue-based methods with real-valued priorities can be very 
efficiently implemented using binary search trees, where the operations of 
insertion, deletion and minimum/maximum location are straightforward In the 
preferred embodiment of the invention, the implementation is related to the 
description in L. Garrido, P. Salembier and L. Garcia, "Extensive Operators in 
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Partition Lattices for Image Sequence Analysis", Signal Processing, (66), 2, 1998, 
pp. 157-180. 

The merging sequence, i.e. a list with the successive merging of 
pairs of video segments, is stored and used to generate a hierarchy. The first level 
5 52 in the hierarchy is defined by key frames from the individual segments 

provided by the video segmentation stage 10. The second level stage 54 in the 
hierarchy is defined by key frames from the clusters generated by the algorithm 
used in the segment merging stage 14. 

For visualization and manipulation, after emptying the hierarchical 
10 queue in the merging algorithm, further merging of video segments is allowed to 
yp build a complete merging sequence that converges into a single segment (i.e., the 

pa i key frame stage 56 represents the whole video clip). The whole video clip 

^ therefore constitutes the third level of the hierarchy. The merging sequence is then 

l r 

Jp represented by a Binary Partition Tree (BPT), which is known to be an efficient 

15 structure for hierarchical representation of visual content In a BPT, each node 

Q (with exception of the leaves, which correspond to the initial shots) has two 

p children. (P. Salembier, L. Garrido, "Binary Partition Tree as an Efficient 

:?? Representation for Filtering, Segmentation, and Information Retrieval", IEEE Intl. 

u 

H 1 Conference on Image Processing, ICIP '98, Chicago, Illinois, October 4-7, 1998.) 

2 0 The BPT also provides the starting point to build a tool for user interaction. 

The tree representation provides an easy-to-use interface for 
visualization and manipulation (verification, correction, augmentation, 
reorganization) of the automatically generated video clusters. Given the 
generality of home video content and the variety of user preferences, manual 
25 feedback mechanisms may improve the generation of video clusters, and 

additionally give users the possibility of actually doing something with their 
videos. 

In a simple interface for displaying the tree representation 50 of the 
merging process, an implementing program would read a merging sequence, and 

3 0 build the binary tree, representing each node of the sequence by a frame extracted 

from each segment. A random frame represents each leaf (shot) of the tree. Each 
parent node is represented by the child random-frame with smaller shot number. 
(Note that the term "random" may be preferred instead of "keyframe" because no 
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effort was done in selecting it). Note also that the representation shown in Figure 
5 is useful to visualize the merging process, identify erroneous clusters, or for 
general display when the number of shots is small, but it can become very deep to 
display when the original number of shots is large. 

A second version of the interface could display only the three 
levels of the hierarchy, i.e., the leaves of the tree, the clusters that were obtained 
as the result of the probabilistic merging algorithm, and the complete- video node. 
This mode of operation should allow for interactive reorganization of the merging 
sequence, so that the user can freely exchange video segments among clusters, 
combine clusters from multiple video clips, etc. Integration of either interface 
with other desired features, like playback of preview sequences when clicking on 
the tree nodes, and VCR capabilities, should be clear to those skilled in this art. 

The invention has been described with reference to a preferred 
embodiment. However, it will be appreciated that variations and modifications 
can be effected by a person of ordinary skill in the art without departing from the 
scope of the invention. Although the preferred embodiment of the invention has 
been described for use with consumer home videos, it should be understood that 
the invention can be easily adapted for other applications, including without 
limitation the summarization and storyboarding of digital movies generally, the 
organization of video materials from news and product-related interviews, health 
imaging applications where motion is involved, and the like. 
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